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ABSTRACT 

' We present the results of a deep, wide-area, optical and near-IR survey of massive high-redshift galaxies. The 

Prime Focus Camera (Suprime-Cam) on the Subaru telescope was used to obtain BRIz' imaging over 2 x 940 
arcmin^ fields, while JKs imaging was provided by the SOFI camera at the New Technology Telescope (NTT) 
for a subset of the area, partly from the ESO Imaging Survey (EIS). In this paper, we report on the properties of 
/T-band-selected galaxies, identified from a total area of ^ 920 arcmin^ to /^vega = 19, of which 320 arcmin^ 
are complete to /Tvega = 20. The BzK selection technique was used to assemble complete samples of about 
^C^' I 500 candidate massive star-forming galaxies (sBzKs) and about 160 candidate massive passively evolving 

. galaxies (pBzKs) at L4 ;$ z ;S 2.5; and the {R — /r)ve2a > 5 color criterion was used to assemble a sample of 

about 850 extremely red objects (EROs). We accurately measure surface densities of 1 .20 ± 0.05 arcmin^^ and 
lO i 0.38±0.03 arcmin^^ for the sBzKs andthepBzKs, respectively. Both sBzKs andpBzKs are strongly clustered, 

■ at a level at least comparable to that of EROs, with pBzKs appearing more clustered than sBzKs. We estimate 

the reddening, star formation rates (SFRs) and stellar masses (M*) for the ensemble of sBzKs, confirming 
p . that to /Tvega 20 typical (median) values are - lO'^ Mq, SFR ~ 190 Mgyr"', and E{B - V) - 0.44. A 

I ' correlation is detected such that the most massive galaxies at z ~ 2 are also the most actively star-forming, an 

O I effect that can be seen as a manifestation of downsizing at early epochs. The space density of massive pBzKs 

at z ~ 1.4 — 2 that we derive is 20%±7% that of similarly massive early-type galaxies at z 0. Adding this 
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^ ■ space density to that of our massive star forming class, sBzKs, in the same redshift range produces a closer 

' comparison with the local early-type galaxy population, naturally implying that we are detecting star formation 

in a sizable fraction of massive galaxies at z > 1 .4, which has been quenched by the present day. Follow-up 
ij optical and near infrared spectroscopy is in progress at the ESO Very Large Telescope (VLT) and at the Subaru 

I telescope, in order to elucidate more thoroughly the formation and evolution of massive galaxies. 

■ Subject headings: galaxies: evolution — galaxies: high-redshift — cosmology: observations — galaxies: 

photometry 



1. introduction 

Despite the recent extraordinary progress in observational 
cosmology and the successful convergence on a single cos- 
mological model, galaxy formation and evolution largely re- 
main an open issue. One critical aspect is how and when 
the present-day most massive galaxies (e.g. elliptical galax- 
ies and bulges with ^ lO'^M©) were built up and what 
type of evolution characterized their growth over cosmic 
time (e.g., Cimatti et al. 2004; Glazebrook et al. 2004, and 
references therein). Indeed, various current renditions of 
the ACDM hierarchical merging paradigm differ enormously 
in this respect, with some models predicting the virtually 
complete disappearance of such galaxies by z = 1 — 2 (e.g., 

' Based on data collected at Subaru Telescope (program S02B-101), which 
is operated by the National Astronomical Observatory of Japan. Also based 
on data collected at the New Technology Telescope (program IDs 169. A- 
0725; 67.A-0244(A); 164.O-056I), which is operated by the European South- 
ern Observatory. 



Cole et al. 2000; Menci et al. 2002; Somerville 2004a) and 
other models predicting a quite mild evolution, more in line 
with observations (e.g., Nagamine et al. 2001; 2005; Granato 
et al. 2004; Somerville et al. 2004b; a direct comparison of 
such models can be found in Fig. 9 of Fontana et al. 2004). 
Moreover, models that provide an acceptable fit to the galaxy 
stellar mass function at z > 1 may differ considerably in the 
actual properties of the galaxies with ^ IO^'M© at z 1, 
with some models predicting very few, if any, passively evolv- 
ing galaxies at these redshifts, at variance with recent findings 
(Cimatti et al. 2004; McCarthy et al. 2004; Daddi et al. 2005a; 
Saracco et al. 2005). 

While various ACDM models may agree with each other 
at z ^ (where they all are tuned) their dramatic divergence 
with increasing redshift gives us powerful leverage to re- 
strict the choice among them, thus aiding understanding of 
the physics of galaxy formation and evolution. Hence, a di- 
rect observational mapping of galaxy evolution through cos- 
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mic time is particularly important and rewarding, especially 
if a significant number of massive galaxies at 1 ^ z 3 can 
be identified and studied. In this regard, the critical ques- 
tions concern the evolution with redshift of the number den- 
sity of massive galaxies and their star formation histories, 
as reflected by their colors and spectral energy distributions 
(SEDs). These questions have just started to be addressed 
by various spectroscopy projects, such as the K20 survey 
(Cimatti et al. 2002a, 52 arcmin^), the Hubble Deep Fields 
(HDFs; Ferguson et al.2000, 5.3 arcmin^ in the HDF-North 
and 4.4 arcmin^ in the HDF-South), the Great Observato- 
ries Origins Deep Survey (GOODS; Giavalisco etal. 2004, 
320 arcmin^ in the North and South fields combined) the 
HST/ACS Ultra Deep Field (S. Beckwith et al. 2006, in prepa- 
ration; 12 arcmin^), the Gemini Deep Deep Survey (Abraham 
et al.2004, 121 arcmin^), and the extension down to z 2 
of the Lyman break galaxy (LBG) project (Steidel et al. 2004, 
^ 100 arcmin^). However, massive galaxies are quite rare 
and likely highly clustered at all redshifts, and hence small 
areas such as those explored so far are subject to large cos- 
mic variance (Daddi et al. 2000; Bell et al. 2004; Somerville 
et al. 2004c). Therefore, although these observation have 
demonstrated that old, passive and massive galaxies do exist 
in the field out to z ~ 2, it remains to be firmly established how 
their number and evolutionary properties evolve with redshift 
up to z 2 and beyond. 

To make a major step forward we are undertaking fairly 
deep, wide-field imaging with the Suprime-Cam on Subaru 
of two fields of 940 arcmin^ each for part of which near- 
IR data are available from ESO New Technology Telescope 
(NTT) observations. The extensive imaging has supported 
the spectroscopic follow-up with the VLT and the Subaru 
telescopes, for which part of the data have already been se- 
cured. The prime aim of this survey is to understand how and 
when the present-day massive galaxies formed, and to this 
end, the imaging observations have been optimized for the 
use of optical/near-IR multi-color selection criteria to identify 
both star-forming and passive galaxies at z ~ 2. 

Color criteria are quite efficient in singling out high redshift 
galaxies. The best-known example is the dropout technique 
for selecting LBGs (Steidel et al. 1996). Besides targeting 
LBGs, color criteria have also been used to search for pas- 
sively evolving galaxies at high redshifts, such as extremely 
red objects (EROs) at redshifts z ^ 1 (Thompson et al. 1999; 
McCarthy 2004) and distant red galaxies (DRGs) at redshifts 
z ^ 2 (Franx et al. 2003). 

Recently, using the highly complete spectroscopic redshift 
database of the K20 survey, Daddi et al. (2004a) introduced 
a new criterion for obtaining virtually complete samples of 
galaxies in the redshift range 1.4 <,z;S2.5, based on B, z 
and Ks^ imaging: star-forming galaxies are identified requir- 
ing BzK = (z - -^^)ab - (B - z)ab > —0.2 (for convenience, 
we use the term sBzKs for galaxies selected in this way); and 
passively evolving galaxies at z ^ 1.4 requiring BzK < —0.2 
and (z — K)ab > 2.5 (hereafter pBzKs). This criterion is red- 
dening independent for star-forming galaxies in the selected 
redshift range, thus allowing us also to select the reddest most 
dust-extinguished galaxies, together with those that are old 
and passively evolving. This should allow for a relatively un- 
biased selection of z 2 galaxies within the magnitude limit 
of the samples studied. 

In this paper observations, data reduction and galaxy pho- 

^ hereafter K band for short 



tometry are described, together with the first results on K- 
band selected samples of distant, high redshift massive galax- 
ies. Compared to optical, the near-IR selection (in particular 
in the K band) offers several advantages, including the rel- 
ative insensitivity of the k-corrections to galaxy type, even 
at high redshift, the less severe dust extinction effects, the 
weaker dependence on the instantaneous star formation ac- 
tivity, and a tighter correlation with the stellar mass of the 
galaxies. Therefore, the study of faint galaxy samples se- 
lected in the near-infrared have long been recognized as ideal 
tools to study the process of mass assembly at high redshift 
(Broadhurst et al. 1992; Kauffmann & Chariot 1998; Cimatti 
et al. 2002a). 

The paper is organized as follows: Section 2 describes the 
observations and the data reduction. Section 3 discusses the 
photometric calibration of the images. Section 4 presents the 
selection and number counts for EROs, sBzKs, and pBzKs. 
Section 5 presents the analysis of the clustering of field galax- 
ies, EROs, sBzKs, and pBzKs. The properties of sBzKs are 
presented in Section 6. Finally, a brief summary is presented 
in Section 7. Throughout the paper, we use the Salpeter IMF 
extending between 0.1 and 100 Mp, and a cosmology with 
l}A=0.7,rjM = 0.3,and/!==//o(kms-i Mpc-i)/100 = 0.71. 
For the sake of comparison with previous works, magnitudes 
and colors in both AB and Vega systems have to be used.^ 

2. OBSERVATIONS 

Two widely separated fields were imaged as a part of our 
survey: one centered at a(J2000)= ll''24"'50\ 5(J2000)= 
-2r42'00" (hereafter Deep3a-F), and the the second, the 
so-called "Daddi field" (hereafter Daddi-F; Daddi et al. 2000) 
centered at a(J2000)= 14''49'"29\ (5(J2000)= 09°00'00". 
Details of the optical and near-IR observations are shown in 
Table 1 . Figure^shows the layout of the two areas observed. 

2.1. Near-IR imaging and data reduction 

Infrared observations in the near-infrared passband J and 
Ks were obtained using the SOFI camera (Moorwood, Cuby 
& Lidman 1998) mounted on the New Technology Telescope 
(NTT) at La Silla. SOFI is equipped with a Rockwell 1024^ 
detector, which, when used together with its large field objec- 
tive, provides images with a pixel scale of 0".29 and a field of 
view of ^ 4.9 x 4.9 arcmin^. 

Deep3a-F is part of the ESO Deep PubHc Survey (DPS) 
carried out by the ESO Imaging Survey (EIS) program (Ren- 
zini & da Costa 1999) (see http://www.eso.0rg/science/eis/1. 
The Deep3a SOFI observations cover a total area of about 
920 arcmin^ in the K band, most at the relatively shallow lim- 
its of /^vega ~ 19.0-19.5. About 320 arcmin^, the region used 
in the present paper, have much deeper integrations with a 
minimum 3600s per sky pixel (and up to 2 hr) reaching to 
^Vega ^ 20 and /vega ^ 22. The NTT J- and /T-band images 
of Deep3a-F were retrieved from the ESO Science Archive 
and reduced using the EIS/MVM pipeline for automated re- 
duction of optical/infrared images (Vandame 2002). The soft- 
ware produces fully reduced images and weight-maps carry- 
ing out bias subtraction, flat-fielding, de-fringing, background 
subtraction, first-order pixel-based image stacking (allowing 
for translation, rotation and stretching of the image) and as- 
trometric calibration. Mosaicking of individual SOFI fields 

^ The relevant conversions between Vega and AB magnitudes for this paper 
are Bab =Bvcga-0.08, Rab =Svcga + 0.22, zab = zvcga + 0.53, and^AB = 

/fVcga + 1-87. 
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Table 1 
Journal of observations. 



Filter 


Telescope 


Obs. date 


Exps.'' 


Seeing 




Area 








(sec) 


(") 


(mag) 


(arcmin^) 


Deep3a-F 


B 


Subai'u 


Mar.5, 03 


3900 


0.77 


27.4 


940 


Rc 


Subaiu 


Mar.4-5, 03 


7320 


0.85 


26.9 


940 


I 


Subaru 


Mar.4-5, 03 


5700 


0.77 


26.5 


940 


z' 


Subaiu 


Mar.4-5, 03 


9900 


0.80 


26.0 


940 


J 


NTT 


Jan.OO-Feb.Ol 


3600 


0.76 


23.4 


320 


Ks 


NTT 


Jan.OO-Feb.Ol 


4800 


0.76 


22.7 


320 


Daddi-F 


B 


Subaru 


Mar 5, 03 


1500 


0.75 


27.0 


940 




WHT 


May 19-21, 98 


3600 


0.70 


25.6 


715 


I 


Subara 


Mar. 5, 03 


1800 


0.90 


26.0 


940 


z' 


Subaru 


Mar 4-5, 03 


2610 


0.80 


25.5 


940 




NTT 


Mar. 27-30, 99 


720 


0.90 


21.5 


600 



"Exposure value for K-band images are "typical values"; see text. 
''The limiting magnitude (in AB) is defined as the brightness corresponding 
to 5 fT on a 2" diameter aperture, 
'^i?— and K-hand data of Daddi-F are described in Daddi et al. (2000). 
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222.200 222.300 222.400 

RA(2000) 
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Fig. 1. — Composite pseudo-color images of Daddi-F (a) and Deep3a-F (b). The RGB colors ai'e assigned to z-, /-, and B-band images, 940 arcmin^, 
respectively. The green area outlined near the center of the images is the field where K-hand images have been obtained by NTT (600 ai'cmin' ai'ea for Daddi-F 
and 320 arcmin^ area for Deep3a-F). 



was based on the astrometric solution. 

Photometric calibration was performed using standard stars 
from Persson et al. (1998), and the calibration performed as 
linear fits in airmass and color index whenever the airmass and 
color coverage allowed for it. The reduced NTT /iT-band data 
and WHT /?-band data for Daddi-F were taken from Daddi 
et al. (2000). The average seeing and the size of the final coad- 
ded images are reported in Table 1 . 

2.2. Optical imaging and data reduction 

Deep optical imaging was obtained with the Prime Fo- 
cus Camera on the Subaru Telescope, Suprime-Cam, which 



with its 10 2k x4k MIT/LL CCDs covers a contiguous area 
of 34' X 27' with a pixel scale of 0."202 pixel"' (Miyazaki 
et al. 2002). Deep3a-F was observed with the four standard 
broad-band filters, B, Rc (hereafter R band), /, and z' (hereafter 
z band) on the two nights of 2003 March 4-5 with 0".7 - 0".9 
seeing. During the same nights also Daddi-F was imaged in B, 
I, and z to a somewhat shallower magnitude limit to match the 
shallower R and K data from Daddi et al. (2000). A relatively 
long unit exposure time of several hundred seconds was used 
in order to reach background-noise-dominated levels. For this 
reason bright stars, which are saturated in the optical images, 
have been excluded by subsequent analysis. 
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During the same nights the photometric standard-star field 
SA95 (Landolt 1992) was observed for -R, and /-band 
flux calibration, and the SDSS standard-star fields SA95- 
190 and SA95-193 were observed for z-band flux calibration 
(Smith et al. 2002). 

The Subaru imaging was reduced using the pipeline pack- 
age SDFRED (Yagi etal.2002; Ouchi etal.2004). The 
package includes overscan correction, bias subtraction, flat- 
fielding, correction for image distortion, PSF matching (by 
Gaussian smoothing), sky subtraction, and mosaicking. Bias 
subtraction and flat fielding were processed in the same man- 
ner as for the conventional single chip CCD. 

In mosaicking, the relative positions (shifts and rotations) 
and relative throughput between frames taken with different 
CCDs and exposures are calculated using stars common to 
adjacent frames and running SExtractor (Bertin & Arnout 
1996) wifli an S/N = 10 flireshold. 

3. PHOTOMETRY 

We obtained K-selected catalogs of objects in our sur- 
vey by detecting sources in the K-band mosaics. For the 
Daddi-F we used the sample of /^-selected galaxies defined 
in Daddi et al. (2000). SExtractor (Bertin & Arnouts 
1996) was used to perform the image analysis and source 
detection in Deep3a-F. The total magnitudes were then de- 
fined as the brightest between the Kron automatic aperture 
magnitudes and the corrected aperture magnitude. Multicolor 
photometry in all the available bands was obtained by run- 
ning SExtractor in double image mode, after aligning all 
imaging to the K-band mosaic. Colors were measured using 
2" diameter aperture magnitudes, corrected for the flux loss 
of stars. The aperture corrections were estimated from the 
difference between the SExtractor Kron automatic aper- 
ture magnitudes (MAG_AUTO) and the 2" aperture magni- 
tudes, resulting in a range of 0.10 - 0.30 mag, depending on 
the seeing. All magnitudes were corrected for Galactic ex- 
tinction (Ab = 0.18 and 0.13 for Deep3a-F and Daddi-F, re- 
spectively) taken from Schlegel et al. (1998), using the empir- 
ical selective extinction function of Cardelli et al. (1989), with 
Rv =Av/E{B-V) ^3.1. 

In Deep3a-F we selected objects to AVega < 20, over a total 
sky area of 320 arcmin^. Simulations of point sources show 
that in all the area the completeness is well above 90% at these 
K-band levels. We recall that objects in Daddi-F were selected 
to completeness limits of AVega < 18.8 over an area of 700 
arcmin^ (but we limit our discussion in this paper to the 600 
arcmin^ covered by the Subaru observations) and to AVega < 
19.2 over a sub-area of 440 arcmin^ (see Daddi et al. 2000 for 
more details). 

The total area surveyed, as discussed in this paper, there- 
fore ranges from a combined area of 920 arcmin (Daddi-F 
and Deep3a-F) at A^/ega < 18.8 to 320 arcmin^ (Deep3a-F) at 
AVega < 20. 

Objects in Deep3a-F and Daddi-F were cross-correlated 
with those available from the 2MASS survey (Cutri 
et al. 2003) in the J and K bands, resulting in good photomet- 
ric agreement at better than the 3% level. In order to further 
verify the photometric zero points we checked the colors of 
stars (Fig. 13, selected from the BzK diagram following Daddi 
et al. (2004a; see Sect l3.1> . From these color-color planes, we 
find that the colors of stellar objects in our data are consis- 
tent with those of Pickles (1998), with offsets, if any, of < 0.1 
mag at most. Similar agreement is found with the Lejeune 



etal.(1997) models. 

Figure 13 shows a comparison of K-band number counts in 
our survey with a compilation of literature counts. No attempt 
was made to correct for different filters (A^ or K). No correc- 
tions for incompleteness were applied to o ur da ta, and we ex- 
cluded the stars, using the method in Sect. 13.11 The filled cir- 
cles and filled squares correspond, respectively, to the counts 
of Deep3a-F and Daddi-F. As shown in the figure, our counts 
are in good agreement with those of previous surveys. 

3.1. Star-galaxy separation 

Stellar objects are isolated with the color criterion (Daddi 
et al. 2004a) (z - K)ab < 0.3(B - z)ab - 0.5. In Fig. H we 
compare the efficiency of such a color-based star-galaxy clas- 
sification with the one based on the SExtractor parameter 
CLASS_STAR, which is based on the shape of the object's 
profile in the imaging data. It is clear that the color classi- 
fication is superior, allowing us to reliably classify stars up 
to the faintest limits in the survey. However, SExtractor 
appear to find resolved in the imaging data a small fraction 
of objects that are color-classified as stars. Most likely, these 
are blue galaxies scattered into the stellar color boundaries by 
photometric uncertainties. 

4. CANDIDATES OF SBzKS, pBzKS AND EROS 

In this section, we select sBzKs, pBzKs and EROs in the 
Deep3a-F and Daddi-F, using the multicolor catalog based 
on the NIR K-band image (see Sect.|3- Forthcoming papers 
will discuss the selection of DRGs (J — K > 2.3 objects) and 
LBGs, using the BRIzJK photometry from our database. 

4. 1 . Selection of sBzKs and pBzKs 

In order to apply the BzK selection criteria consistently 
with Daddi et al. (2004a), we first accounted for the differ- 
ent shapes of the filters used, and applied a correction term 
to the B-band. The B-band filter used at the Subaru tele- 
scope is significantly redder than that used at the VLT by 
Daddi et al. (2004a). We then carefully compared the stellar 
sequence in our survey to that of Daddi et al. (2004a), using 
the Pickles (1998) stellar spectra and the Lejeune et al. (1997) 
corrected models as a guide, and applied small color terms to 
B — z andz — (smaller than ^0.1 mags in all cases), in order 
to obtain a fully consistent match. In the following we refer to 
BzK photometry for the system defined in this way, consistent 
with the original BzK definition by Daddi et al. (2004a). 

Figure |5l shows the BzK color diagram of K-selected ob- 
jects in Deep3a-F and Daddi-F. Using the color criterion from 
Daddi et al. (2004a), BzK = (z - K)pji - (B - z)ab > -0.2, 
387 galaxies with /Tvega < 20 were selected in Deep3a-F as 
sBzKs, which occupy a narrow range on the left of the solid 
line in Fig.|51i. Using BzK < -0.2 and (z-/:)ab > 2.5, 121 
objects were selected as candidate pBzKs, which lie in the 
top-right part of Fig.|5t. To /Tvega < 20 a surface density of 
1 .20 ± 0.05 arcmin"^ 9 33 _^ o.03 arcmin"^ derived 

separately for sBzKs and pBzKs, respectively (Poisson's er- 
rors only). The surface density of sBzKs is larger but still 
consistent within la with the 0.91 ±0.13 arcmin^^ found in 
the 52 arcmin^ of the K20 field (Daddi et al. 2004a), and with 
the 1.10 ±0.08 arcmin"2 f^^^^ (^e GOODS North field 
(Daddi et al. 2005b). Instead, the surface density of pBzKs 
recovered here is significantly larger than what found in both 
fields. This may well be the result of cosmic variance, given 
the strong clustering of pBzKs (see Section|5j, and their lower 
overall surface density. 
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Fig. 2. — Optical two-color plots of A'-selected objects in our survey: the left panels show B — R vs. R — I colors; the right panels ha ve B — I vs. R — z colors 
(all in AB scale). Galaxies are shown as filled points and stars with green asteiisks (based on the BzK color star-galaxy seperation, Sect. l3.ll . [See the electronic 
edition of the Journal for the color version of this figure.} 



Using the same criteria, we select sBzKs and pBzKs in the 
Daddi-F field. In Daddi-F 108 sBzKs and 48 pBzKs are se- 
lected, and they are plotted in Fig.|31;. The density of sBzKs 
and pBzKs in Daddi-F is consistent with that in Deep3a, if 
limited at A'vega ^19. 

4.2. Selection ofEROs 

EROs were first identified in /T-band surveys by Elston, 
Rieke, & Rieke (1988), and are defined here as objects hav- 
ing red optical-to-infrared colors such that {R — ^r)vega > 
5-6, corresponding to {R — K)p^z > 3.35^.35. EROs are 
known to be a mixture of mainly two different populations 
at z ^ 0.8; passively evolving old elliptical galaxies and dusty 
starburst (or edge-on spiral) galaxies whose UV luminosities 
are strongly absorbed by internal dust (Cimatti et al. 2002b; 
Yan et al. 2004, Moustakas et al. 2004). In Daddi-F EROs 
were selected and studied by Daddi et al. (2000) using various 
R — K thresholds. In order to apply a consistent ERO selection 
in Deep3a-F, we considered the filter shapes and transmission 
curves. While the same K-band filters (and the same tele- 
scope and instrument) were used for K-band imaging in the 
two fields, the R-band filters used in the two fields differ sub- 
stantially. In the Daddi-F the WHT R-band filter was used, 
which is very similar, e.g., to the R-band filter of FORS at 
the VLT used by the K20 survey (Cimatti et al. 2002a). The 



SubaruH-Suprime-Cam /?-band filter is much narrower than 
the above, although it does have a very close effective wave- 
length. As a result, distant z ^ I early-type galaxy spectra 
as well as M-type stars appear to have much redder R — K 
color, by about «0.3 mag, depending on exact redshift and 
spectral shape. Therefore, we selected EROs in Deep3a-F 
with the criterion Rsuburu — K > 3.7 (AB magnitudes), cor- 
responding closely to Rwht — K > 3.35 (AB magnitudes) or 
RwHT —K>5 (Vega magnitudes). In Deep3a-F, 513 EROs 
were selected to /Tvega < 20, and they are plotted in Fig. [Sj) 
with solid red triangles, for a surface density of 1 .6 arcmin^^. 
To the same AVega < 20 limit, this agrees well with the density 
found, e.g., in the K20 survey (~ 1 .5 arcmin^^), or in the 180 
arcmin^ survey by Georgakakis et al. (2005). In the Daddi-F, 
337 EROs were selected with the criterion Rwht — K > 3.35, 
consistent with what done in Deep3a, and are plotted in 
Fig. I^Jl as red solid triangles. The surface density of EROs 
in both fields at /Tvega < 1 8 .4 can be compared, with overall 
good consistency, to the one derived from the large 1 deg^ 
survey by Brown et al. (2005). 

The peak of the EROs redshift distribution is at z 1 (e.g., 
Cimatti et al. 2002a). By looking at the BzK properties of 
EROs we can estimate how many of them lie in the high-z tail 
z > 1 .4, thus testing the shape of their redshift distribution. 
In the Deep3a-F to A'vega < 20 some 90 of the EROs are also 
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Fig. 3. — Differential K band galaxy counts from Deep3a-F and Daddi-F, compared with a compilation of results taken from various sources. [See the 
electronic edition of the Journal for the color version of this figure.} 



sBzKs, thus likely belong to the category of dusty starburst 
EROs at z > 1.4, while 121 EROs are classified as pBzKs. 
Totally, ^41% of EROs are selected with the BzK criteria, 
thus expected to lie in the high-z tail (z > 1 .4) for ^Tve^a < 20 
sample. This result is consistent with the value of 35% found 
in the 52 arcmin^ of the K20 field (Daddi et al. 2004a), and 
with the similar estimates of Moustakas et al. (2004) for the 
full GOODS-South area. In the Daddi-F, to /Tvega < 19.2, 49 
of the EROs are also sBzKs; and 48 of them are also pBzKs. 
About 29% of EROs at iiTvega < 19.2 are in the high-z tail at 
z> 1.4. 

4.3. Number counts of EROs, sBzKs, and pBzKs 

Simple surface densities provide limited insight into the 
nature of different kind of galaxies. However, number- 
magnitude relations, commonly called number counts, pro- 
vide a statistical probe of both the space distribution of galax- 
ies and its evolution. For this reason, we derived K-bwd 
differential number counts for EROs, sBzKs and pBzKs in 
our fields, and plotted them in Figure |6l The differential 
number counts in 0.5 mag bins are shown in Table |2l Also 
shown are the number counts for all /T-selected field galax- 
ies in Deep3a-F (circles with solid line) and Daddi-F (squares 
with dot-dashed line), with the same as in Fig.l^for compar- 
ison. A distinguishable characteristic in Fig.|6lis that all of 
high redshift galaxies (EROs, sBzKs, and pBzKs) have faint 



NIR apparent magnitudes (.fiTvega ^17 mag), and the slopes of 
the counts for EROs, sBzKs and pBzKs are steeper than that 
of the full /T-selected sample. 

The open squares with a solid line in Fig.|6lshows the num- 
ber counts for sBzKs in Deep3a-F. The fraction of sBzKs in 
Deep3a-F increases very steeply towards fainter magnitudes. 
The triangles with a solid line and crosses with a solid line 
show, respectively, the number counts of EROs and pBzKs in 
Deep3a-F. The open squares with a dot-dashed line in Fig.|6l 
show the number counts for sBzKs in Daddi-F. The counts of 
sBzKs in Daddi-F are almost identical to those in Deep3a-F, 
to their limit of &/ega ~ 19. The triangles with a dot-dashed 
line and crosses with a dot-dashed line in Fig.|6lshow respec- 
tively the number counts for EROs and pBzKs in Daddi-F. 

For EROs, the slope of the number counts is variable, be- 
ing steeper at bright magnitudes and flattening out toward 
faint magnitudes. A break in the counts is present at /Tvega ~ 
18.0, very similar to the break in the ERO number counts 
observed by McCarthy et al. (2001) and Smith et al. (2002). 
The pBzKs number counts have a similar shape, but the 
break in the counts slope is apparently shifted ^ 1-1.5 mag 
fainter There are indications that EROs and pBzKs have 
fairly narrow redshift distributions: peaked at z 1 for EROs 
(Cimatti et al. 2002b; Yan et al. 2004; Doherty et al. 2005) and 
at z « 1 .7 for pBzKs (Daddi et al. 2004b; 2005a). The number 
counts might therefore be direct probes of their respective lu- 
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Fig. 4. — Star/galaxy separation. Left: K-band magnitude vs. the stellarity index parameter (CLASS_STAR) from SExtractor for objects (small 
dotted points) in Deep3a-F. The dashed lines are CLASS_STAR = 0.95 and A'vega = 18.0; objects with {z, — A')ab — 0.3(B — z)ab < —0.5 are plotted as open 
circles. Right: B — z against z — K for objects in Deep3a-F. objects with CLASS_STAR > 0.95 and ^Tvcga < 18.0 are plotted as triangles, and those with 
(z — X')ab — 0.3{B — z)ab < —0.5 are plotted as open circles. Stars con'espond to stellar objects given in Pickles (1998). The dot-dashed line, {z — K) = 
0.3 (B — z) — 0.5 , denotes the boundary between stars and galaxies adopted in this study; an object is regarded as a star, if it is located below this line. [See the 
electronic edition of the Journal for the color version of this figure.} 
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Table 2 

Differential number counts in 0.5 magnitude bins of EROs, sBzKs, and pBzKs in Deep3a-F and Daddi-F. 





Deep3a-F in log (N/deg-/0.5mag) 








Daddi-F in 


log (N/deg-/0.5mag) 






K bin center 


Galaxies 


EROs 


sBzKs 


pBzKs 


K bin center 


Galaxies 


EROs 


sBzKs 


pBzKs 


16.75 


2.981 


1.353 






16.75 


2.888 


1.254 






17.00 


3.109 


1.654 






17.00 


3.023 


1.555 






17.25 


3.213 


1.830 






17.25 


3.165 


1.891 






17.50 


3.337 


2.198 






17.50 


3.323 


2.120 






17.75 


3.470 


2.483 


1.353 




17.75 


3.426 


2.321 


1.254 




18.00 


3.565 


2.675 


1.654 


1.052 


18.00 


3.498 


2.590 


1.622 


1.078 


18.25 


3.678 


2.822 


2.006 


1.830 


18.25 


3.569 


2.786 


1.923 


1.379 


18.50 


3.764 


3.025 


2.228 


2.353 


18.50 


3.669 


2.888 


2.209 


1.891 


18.75 


3.802 


3.138 


2.529 


2.467 


18.65 


3.708 


2.990 


2.342 


2.175 


19.00 


3.859 


3.145 


2.724 


2.596 


19.00 


3.783 


3.063 


2.834 


2.461 


19.25 


3.911 


3.162 


2.971 


2.608 












19.50 


3.983 


3.201 


3.228 


2.675 












19.75 


4.072 


3.297 


3.470 


2.759 













minosity function, and the shift in the counts is indeed consis- 
tent with the different typical redshifts of the two population 
of galaxies. 

The counts of sBzKs have roughly the same slope at all K- 
band magnitudes. This is consistent with the much wider red- 
shift distribution of this class of galaxies. However, we expect 
that at bright magnitudes AGN contamination might be more 
relevant than at Ky^^^ ~ 20. Correcting for this, the counts 
of non-AGN sBzKs galaxies at bright magnitudes might be 
intrinsically steeper 

5. CLUSTERING OF K-SELECTED GALAXY POPULATIONS 

Measuring the clustering of galaxies provides an additional 
tool for studying the evolution of galaxies and the formation 



of structures. In this section we estimate over the two fields 
the angular correlation of the general galaxy population as 
well as of the EROs, sBzKs and pBzKs. 

In order to measure the angular correlation function of the 
various galaxy samples, we apply the Landy & Szalay tech- 
nique (Landy & Szalay 1993; Kerscher et al. 2000), follow- 
ing the approach already described in Daddi et al. (2000), to 
which we refer for formulae and definitions. 

5.1. Clustering of the K-selected field galaxies 

In our analysis a fixed slope (5 = 0.8 was assumed for the 
two-point correlation function [w{6) =Ax 10^*^]. This is 
consistent with the typical slopes measured in both faint and 
bright surveys and furthermore makes it possible to directly 



8 



X. Kong et al. 





0246 0246 
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Fig. 5. — Two-color (z — K)xb vs {B — z)ab diagram for the galaxies in the Deep3a-F and Daddi-F fields. Galaxies at high redshifts are highlighted. The 
diagonal solid line defines the region BzK = (z — A^)ab ^ — z)ab > —0.2 that is efficient to isolate z > 1.4 star forming galaxies (sBzKs). The horizontal 
dot-dashed line further defines the region (z — X^) ab > 2.5 that contains old galaxies at z > 1 .4 (pBzKs). The dashed fines seperate regions occupied by stars and 
galaxies. Filled stars show objects classified as stars, having (z — ^^^)ab — 0.3 (B — z)ab < —0.5; open stars show stellar objects from the K20 survey (Daddi et al. 
2005a); squares represent sBzKs; circles represent pBzKs; triangles represent galaxies with {R — K)xb > 3.35 (EROs). Galaxies lying out of the BzK regions, 
thus likely having in general redshifts less than 1.4, are simply plotted as black points. 



compare our results with the pubhshed ones that were ob- 
tained adopting the same slope (see Daddi et al. 2000 for more 
details). 

In Figures and |8] the bias-corrected two-point correla- 
tion functions w{9) of Daddi-F and Deep3a-F are shown as 
squares; the bins have a constant logarithmic width (A log 6 = 
0.1), with the bin centers ranging from 5.6" to 11. 8' for both 
fields. These values of 9 are large enough to avoid prob- 
lems of under-counting caused by the crowding of galaxy 
isophotes and yet are much smaller than the extent of the in- 
dividual fields. The dashed line shows the power-law cor- 
relation function given by a least square fit to the measured 
correlations. We clearly detect a positive correlation signal 
for both fields with an angular dependence broadly consis- 
tent with the adopted slope S — 0.8. The derived cluster- 
ing amplitudes (where A is the amplitude of the true angu- 
lar correlation at 1°) are presented in the third column of Ta- 
ble|3]and of Table|4]for Daddi-F and Deep3a-F, respectively, 
and shown in Fig. |9l We compare our result on Daddi-F 
with those previously reported by Daddi et al. (2000), since 
we use the same /T-band images for this field. The present 
^T-selected galaxy sample is however different from that of 
Daddi et al. (2000) in thta: (1) the area (600 arcmirr-) is 
smaller than that (701 arcmin^) in Daddi et al. (2000), because 



we limit ours to the area with all the BRIzK-hand data; and (2) 
the star-galaxy separation was done by different methods. In 
Daddi et al. (2000) the SExtractor CLASS_STAR mor- 
phological parameter was used, while the photometric crite- 
rion of Daddi et al. (2004a) is used here. Column lOof Table|3] 
lists the clustering amplitudes of all the TiT-selected galaxies in 
Daddi et al. (2000). For /Tvega ==18.5 and fainter bins the A- 
values in the two samples are in fair agreement, but for the 
brightest bin (/Tvega = 18) the A-value in Daddi et al. (2000) is 
1.3 X 10"^ 45% smaller than that found here (2.36 x 10"^), 
probably mainly due to t he m ore efficient star-galaxy separa- 
tion employed here (Sect. im . Apart from this small discrep- 
ancy, we find good agreement between our results and those 
of Daddi et al. (2000) and Oliver et al. (2004). 

The clustering amplitudes in Deep3a-F tend to be slightly 
but systematically smaller than in Daddi-F, which is likely 
caused by the intrinsic variance among different fields, de- 
pending on the survey geometry, surface density and cluster- 
ing properties (see Sect. 15. 5> . 

We find a smooth decline in amplitude with K-band 
magnitude that is consistent with the results from Roche 
et al. (2003), Fang et al. (2004), and Oliver et al. (2004). The 
decline is not as steep as in the range 15 < Ky^ga < 18 (see 
e.g. Roche et al. 1999). The flattening, which extends up 
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Fig. 6. — K-band differential galaxy number counts for EROs, sBzKs and pBzKs, compared with the ^T-selected field galaxies as in Fig.|3] The sohd curves 
show the number counts for objects in Deep3a-F, and the dot-dashed curves show the number counts for objects in Daddi-F. Triangles, open squares and crosses 
show the number counts for EROs, sBzKs and pBzKs, respectively. The circles and squares show the iT-selected field galaxies in I3eep3a-F and Daddi-F, 
respectively. 



Table 3 

Clustering amplitudes for the K-selected sample, EROs, sBzKs, and pBzKs in Daddi-F. 



K limit' 




Galaxies 




ERO.S 




sBzKs 




pBzKs 


Daddi et al. (2000) 


num. 


ALIO--']" 


num. 


A[10--'] 


num. 


A[10--] 


num. 


A[10--'] 


Gal.= 


EROs' 


18.0 


978 


2.36±0.94 


51 


23.60±4.18 


7 




1 




1.3 


24 


18.5 


1589 


2.16±0.40 


132 


22.00±2.82 


21 


24.00±9.80 


5 




1.6 


22 


18.8 


2089 


1.91±0.30 


228 


14.60±1.64 


43 


18.90±7.52 


20 


24.70±9.92 


1.5 


14 


19.2 


2081 


1.68±0.28 


264 


13.20±1.26 


92 


11.50±6.61 


40 


22.90±7.63 


1.6 


13 



"The area for Kvega < 18.8 mag is about 600 arcmin^, and for Kvega = 19.2 
it is about 440 arcmin^ . 



''The amplitude of the true angular correlation at 1°, the value of the C is 
5.46 and 5.74 for the whole and the deeper area, respectively. 

'^The last two columns show the clustering amplitudes for the K-selected 
galaxies and the EROs in Daddi et al. (2000). 



Table 4 

Clustering amplitudes for the K-selected sample, EROs, sBzKs, and pBzKs in Deep3a-F. 



K limit 




Galaxies 




EROs 




sBzKs 




pBzKs 


num. 


A[10--']" 


num. 


A[10-'] 


num. 


A[10-^] 


num. 


A[10-^] 


18.5 


997 


1.96±1.04 


95 


14.70±2.20 


13 




8 




18.8 


1332 


1.65±0.76 


166 


9.29±1.60 


27 


29.50±5.80 


26 


40.90±7.55 


19.5 


2284 


1.24±0.41 


340 


4.89±0.78 


129 


6.70±3.14 


71 


21.40±3.00 


20.0 


3333 


1.14±0.28 


513 


4.25±0.52 


387 


4.95±1.69 


121 


10.40±2.83 



"The amplitude of the true angular coiTelation at 1°, the value of the C is 
6.63. 
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to K ^ 22-24 (Carlberg et al. 1997; Daddi et al. 2003) has 
been interpreted as due to the existence of strongly clustered 
K-selected galaxy populations extending to redshifts z 1- 
3 (Daddi et al. 2003). Beside the well known z « 1 EROs 
strongly clustered galaxy populations, other populations with 
high angular clustering indeed exist with redshift extending to 
z ~ 2.5 at least, as discussed in the next sections. 

5.2. Clustering of the EROs 

We estimate the clustering properties of the EROs, using 
the large sample of EROs derived from our two fields. Fig- 
uresllOh) andllOb) clearly show, from both fields, that the sky 
distribution of EROs is very inhomogeneous, as first noted by 
Daddi et al. (2000). The two-point correlation functions are 
shown in Fig.0and|8]for Daddi-F and Deep3a-F, and the dot- 
dashed lines show the power-law correlation function given 
by a least squares fit to the measured values. The correlations 
are well fitted by a (5 = 0.8 power law. 

A strong clustering of the EROs is indeed present at all 
scales that could be studied, and its amplitude is about one or- 
der of magnitude higher than that of the field population at the 
same /Tvega limits, in agreement with previous findings (Daddi 
et al. 2000; Firth et al. 2002; Brown et al. 2005; Georgakakis 
et al. 2005). The derived clustering amplitudes are reported in 
column (5) of Tableland Table|4]for Daddi-F and Deep3a-F, 
respectively. The amplitudes shown in Fig.|5]suggest a trend 
of decreasing strength of the clustering for fainter EROs in 
both fields. 

For the ^^vega < 18.5 and 18.8 mag subsamples of Daddi-F, 
the correlation function signal is significant at the 7 a level 
with clustering amplitudes A = 22.0 x 10^^ and 14.6 x 10^^, 
respectively. For the subsample with /Tvega < 19.2 mag, the 
detected signal is significant at the 10 cr confidence level 
(A = 13.2 X 10"^). In column 11 of Table^we also list the 
clustering amplitudes of EROs in Daddi et al. (2000), which 
are in good agreement with the present findings. Using a 
~ 180 arcmin^ Ks-band survey of a region within the Phoenix 
Deep Survey, Georgakakis et al. (2005) have analyzed a sam- 
ple of 100 EROs brighter than /Tvega = 19 mag, and estimated 
an amplitude A = 1 1 .7 x 10^^, consistent with our results. 

As for the TiT-selected galaxies, the clustering amplitude of 
EROs in Deep3a-F is slightly smaller than that in Daddi-F. 
The clustering amplitudes are A = 9.29 X 10-^ for /iTvega < 
18.8 EROs, andA = 4.25 x 10"^ for/Tvega < 20. EROs, which 
is weaker than the A = 8.7 x 10"^ for /Tvega < 20 EROs in 
Georgakakis et al. (2005). Field-to-field variation is one of 
possi ble reasons for this discrepancy, which is discussed in 
Sect. 1531 

We notice that our results are solid against possible contam- 
ination by stars. Contamination by unclustered populations 
(i.e. stars) would reduce the amplitude of the angular corre- 
lation function by (1 — /)^, where / is the fractional contam- 
ination of the sample. The prime candidates for contamina- 
tion among EROs are red foreground Galactic stars (note that 
stars are not contaminating BzK samples), which can have 
red R K colors. However, we have rejected stars among 
EROs using th e pho tometric criterion for star-galaxy separa- 
tions (see Sect. 13. 1> . Therefore, our clustering measurements 
refer to extragalactic EROs only. 

5.3. Clustering of the star-forming BzKs 

Figures llOt ) and HOB ) display the sky distribution of the 
sBzKs in Daddi-F and Deep3a-F, and shows that also these 



galaxies have a quite inhomogeneous distribution. This is not 
an artifact of variations of the detection limits over the fields, 
because the Monte Carlo simulations show that differences 
in detection completeness within small (4'.0 x 4'.0) areas in 
the image are very small, and that the detection completeness 
does not correlate with the distribution of the sBzKs. 

The resulting angular correlation functions for the sBzKs 
are shown in Fig. and Fig.|8lfor the two fields. Again, a 
slope 5 = 0.8 provides a good fit to the data. The best fit 
values of A are reported in column 7 of Tableland Table|3 

The w{9 = 1°) amplitudes of sBzKs in Daddi-F ai-e 24.0 x 
10-^ 18.9 X 10-^ and 11.5 x 10"^ at^Tvega = 18.5, 18.8 and 
19.2, respectively. For sBzKs in Deep3a-F, the w{6) ampli- 
tudes become 29.5 x 10"^ 6.70 x 10"^ and 4.95 x 10"^ at 
^Vega — 18.8, 19.5 and 20.0, respectively. The sBzKs ap- 
pear to be strongly clustered in both fields, and the clustering 
strength increases with the /T-band flux. Actually, they are 
as strongly clustered as the EROs. Strong clustering of the 
sBzKs was also inferred by Daddi et al. (2004b), by detecting 
significant redshift spikes in a sample of just nine /Tvega < 20 
sBzKs with spectroscopic redshifts in the range 1.7 < z < 
2.3. Albeit without spectroscopic redshifts, the present re- 
sult is instead based on a sample of 500 sBzKs. Adelberger 
et al. (2005) also found that UV-selected, star forming galax- 
ies with /Tvega < 20.5 in the redshift range 1.8 < z < 2.6 are 
strongly clustered. Our results are in good agreement with 
these previous findings. 

5.4. Clustering of the passive BzKs 

Figures 1 10b ) and IIOF ) display the sky distribution of the 
pBzKs in Daddi-F and Deep3a-F, and show that also these 
galaxies have a very inhomogeneous distribution. We then 
derive the angular two-point correlation function of pBzKs, 
using the same method as in the previous subsections. 

The resulting angular correlation functions for the pBzKs 
are shown in Fig. and Fig.|8lfor the two fields. Again, a 
slope S = 0.8 provides a good fit to the data. The best fit 
values of A are reported in column (9) of Tableland Tabled 

The pBzKs appear to be the most strongly clustered galaxy 
population in both fields (with A = 22.9 x 10^^ for /Tvega < 
19.2pBzKs in Daddi-F and A = 10.4 x 10"^^ for /Tvega < 20.0 
pBzKs in Deep3a-F), and the clustering strength increases 
with increasing /T-band flux. 

5.5. K-band dependence and field-to-field variations of 
clustering measurements 

Fig-Elsummarizes the clustering measurements for the pop- 
ulations examined (field galaxies, EROs, and BzKs), as a 
function of the A'-band limiting magnitudes of the samples. 
Clear trends with K are present for all samples, showing that 
fainter galaxies have likely lower intrinsic (real space) clus- 
tering, consistent with the fact that objects with fainter K are 
less massive, or have wider redshift distributions, or both. For 
the EROs (see Sect. 14. 2> we have found evidence that faint 
Kveoa < 20 samples have indeed a higher proportion of galax- 
ies in the z > 1.4 tail, with respect to brighter /Tvega < 19 
objects, and thus a wider redshift distribution. As already 
noted, all color-selected high-redshift populations are sub- 
stantially more clustered than field galaxies, at all the magni- 
tudes probed here. The reason for the stronger angular cluster- 
ing of pBzKs, compared, e.g, to sBzKs or EROs, is likely (at 
least in part) their narrower redshift distributions 1 .4 < z < 2 
(Daddi et al. 2004a; 2005a). In future papers, we will use the 
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Fig. 7. — Observed, bias-corrected two-point coiTelations for the Daddi-F sample of field galaxies (squares), EROs (triangles), sBzKs (filled circles), and 
pBzKs (stars). The error bars on the direct estimator values are Icr errors. To make this plot clear, we show the error bars of EROs and pBzKs only, the error 
bars of field galaxies are smaller, and the error bars of sBzKs are larger than those of EROs. The lines (dashed, dot-dashed, dotted, dash-dot-dotted) show the 
power-law fitted to the w(8). Because of the small number of objects included, some bins were not populated. 



Limber equation, with knowledge of the redshift distributions, 
to compare the real space correlation length of the different 
populations. 

For strongly clustered populations, with angular clustering 
amplitudes A « 10^^, a large cosmic variance of clustering 
is expected, which is relevant also on fields of the size of 
the ones studied here (see Daddi etal. 2001; 2003). There 
are in fact variations in our clustering measurements for high- 



redshift objects between the two fields, sometimes larger than 
expected on the basis of the errors on each clustering mea- 
surements. We remind the reader that realistic external er- 
rors on the angular clustering of EROs, as well as sBzKs 
and pBzKs are likely larger than the Poissonian ones that we 
quote. Following the recipes by Daddi et al. (2001) we esti- 
mate that the typical total accuracy for our measurements of 
A K, 0.02, when including external variance, is on the order of 
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Fig. 8. — Same as FiglTl but for Deep3a-F. 
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5.6. Cosmic variance in the number counts 



The presence of strong clustering will also produce substan- 
tial field to field variations in the number counts. Given the 
available measurements of angular clustering, presented in the 
previous sections, we are able to quantify the expected vari- 
ance in the galaxies counts, following e.g. eq. (8) of Daddi 
et al. (2000). We estimate that, for the Deep3a-F limit of 
^Vega < 20, the integrated numbers of objects are measured 



to ^ 20% precision for EROs and sBzKs, and to 30% preci- 
sion for pBzKs. 

6. PROPERTIES OF SBzKS 

The accurate analysis of the physical properties of high 
redshift galaxies (such as SFR, stellar mass, etc.) requires 
the knowledge of their spectroscopic redshift. VLT VIMOS 
spectra for objects culled by the present sample of sBzKs and 
EROs have been recently secured and are now being analyzed, 
and will be used in future publications. In the meantime, es- 
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timates of these quantities can be derived on the basis of the 
present photometric data, following the recipes calibrated in 
Daddi et al. (2004a), to which we refer for definitions and a 
more detailed discussion of the recipes. 

While errors by a factor of 2 or more may affect individual 
estimates done in this way, when applied to whole popula- 
tion of Bz/T-selected galaxies these estimates should be rela- 
tively robust, on average, because the Daddi et al.'s (2004a) 
relations were derived from a sample of galaxies with spec- 
troscopic redshifts. The estimates presented here thus rep- 
resent a significant improvement on the similar ones provided 
by Daddi et al. (2004a) because of the 6 to 20 times larger area 
probed (depending on magnitude) with respect to the K20 sur- 
vey, which should help to significantly reduce the impact of 
cosmic variance. 

6.1. Reddening and star formation rates 

Following Daddi et al. (2004a), estimates of the reddening 
E{B — V) and SFR for sBzKs can be obtained from the BzK 
colors and fluxes alone (see Daddi et al. 2004a for more de- 
tails). The reddening can be estimated by the B — z color, pro- 
viding a measure of the UV slope. The Daddi et al. 's recipe is 
consistent with the recipes by Meurer et al. (1999) and Kong 
et al. (2004) for a Calzetti et al. (2000) extinction law, based 
on the UV continuum slope. Daddi et al. (2004a) showed 



that for BzK galaxies this method provides E{B — V) with an 
rms dispersion of the residuals of about 0.06, if compared to 
values derived with knowledge of the redshifts, and with the 
use of the full multicolor SED. With knowledge of reddening, 
the reddening corrected B-band flux is used to estimate the 
1500A rest-frame luminosity, assuming an average redshift of 
1.9, which can be translated into SFR on the basis, e.g., of 
the Bruzual & Chariot (2003) models. Daddi et al. (2005b) 
showed that SFRs derived in this way are consistent with ra- 
dio and far-lR based estimates, for the average sBzKs. 

The — y ) and SFR histograms of the sBzKs in Daddi-F 
and Deep3a-F are shown in Fig. ^2 About 95% of the sBzKs 
in Daddi-F (/Tvega < 19.2) have SFR> VOMoyr"', and the 
median SFR is about 370 Mgyr^^ About 90% of the sBzKs 
in Deep3a-F (TiTvega < 20.0) have SFR>7OM0yr-', and the 
median SFR is - 190 Mgyr^^ 

The median reddening for the Ky^o^ < 20 sBzKs is es- 
timated to be E{B — V) = 0.44, consistently with Daddi 
et al. (2004a; 2005b). Of sBzKs 55% have E{B - V) > 0.4, 
the limit at which we estimate the UV-based criteria of Stei- 
del et al. (2004) would fail at selecting z^2 galaxies. There- 
fore, we estimate that ^55% of the z ^ 2 galaxies would 
be missed by the UV criteria. This is similar to, but higher, 
than the 40% estimated by Reddy et al. (2005). The probable 
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reason for the small discrepancy is that Reddy et al. (2005) 
excluded from their sample the reddest sBzKs for which the 
optical magnitudes could not be accurately measured in their 
data. 

6.2. Stellar masses of sBzKs and pBzKs 

Using BC03 models, spectroscopic redshifts of individual 
K20 galaxies, and their whole UBVRIzJHK SED, Fontana 
et al. (2004) have estimated the stellar-mass content for all the 
K20 galaxies (using a Salpeter IMF from 0.1 to 100 Mq, as 
adopted in this paper). The individual mass estimates for 31 
non-AGN sBzKs and pBzKs objects with z > 1.4 have been 
used by Daddi et al. (2004a) to calibrate an empirical relation 
giving the stellar mass for both sBzKs and pBzKs as a func- 
tion of their observed A'-band total magnitude and z~K color. 
This relation allows one to estimate the stellar mass with un- 
certainties on single objects of about 60% compared to the 
estimates based on knowledge of redshifts and using the full 
multicolor SEDs. The relatively small variance is introduced 
by intrinsic differences in the luminosity distance, in the M/L 
ratio for given magnitudes and/or colors. 

The histograms for the stellar mass of the sBzKs derived in 
this way are shown in Fig. II lb and Fig. II IF. About 95% of the 
sBzKs in Daddi-F have > 10^ 'M©, and the median stellar 
mass is 2.0 x lO'^Mg; in Deep3a-F - 40% of the sBzKs have 
M* > lO'^M©, the median stellar mass is ^ 8.7 x lO^^Mo, a 
difference due to the different limiting K magnitude in the two 
fields. 

Using the same method, we also estimate the stellar mass 
of the pBzKs in both fields, and plot them in Fig. Illfc and 
Fig. II IF as the dotted line. The median stellar mass of pBzKs 
in Daddi-F is ~ 2.5 x lO'^M©, and in Deep3a-F is 1.6 x 
1O"M0, respectively. Again, the higher masses for sBzKs 
in Daddi-F compared to Deep3a-F result from the shallower 
/T-band limit. 

It is worth noting that in the /Tvega < 20 sample there 
are barely any pBzKs less massive than 7 x lO'^Mg, (see 
Fig- II IF ), while over 50% of sBzKs are less massive than 
this limit. This is primarily a result of pBzKs being by con- 
struction redder than (z — K)j^q = 2.5, and hence Eq. (6) 
and Eq. (7) in Daddi et al. (2004a) with /Tvega < 20 implies 
M* > 7 x IO'^Mq. Note that above lO^M© (above which 
our sample should be reasonably complete) the numbers of 
sBzKs and pBzKs are similar. We return to this point in the 
last section. 

6.3. Correlation between physical quantities 

Figures and 1 12b) show the correlation between color 
excess E{B — V) and SFR for the sBzKs in Daddi-F and 
Deep3a-F, respectively. The Spearman rank correlation co- 
efficients are r^ = 0.40 for Daddi-F and r^ = 0.57 for Deep3a- 
F. This implies that the SFR is significantly correlated to 
E{B — V), at a > 5ct level, and the reddest galaxies have the 
highest SFR. Part of this correlation can arise from simple 
error propagation, as an overestimate (underestimate) of the 
reddening automatically leads to an overestimate (underesti- 
mate) of the SFR of a galaxy. In Fig.ll2h) and Fig. 112b) the 
arrow shows the resulting slope of the correlated reddening 
errors, which indeed is parallel to the apparent correlation. 
However, the scatter in the original SFR-£'(B — V ) correlation 
[SE{B - y) - 0.06; Daddi et al. 2004a] is much smaller than 
what needed to produce the full correlation seen in Fig.ll2h) 
and 1 12b). We conclude that there is evidence for an intrinsic 



correlation between SFR and reddening for z^2 star-forming 
galaxies, with galaxies with higher star formation rates having 
more dust obscuration. A positive correlation between SFR 
and reddening also exists in the local universe (see Fig. 1 of 
Calzetti 2004), and was also found by Adelberger & Steidel 
(2000) for z ^3 LBGs, on a smaller range of reddening. 

In Figures[T2b) andll2tl). we plot the relation between color 
excess E{B — V) and stellar mass of the sBzKs in Daddi-F 
and Deep3a-F. The Spearman rank correlation coefficient is 
rs — 0.53 for Daddi-F and r^ = 0.63 for Deep3a-F, indicat- 
ing that the correlation between E{B — V) and stellar mass is 
significant at the > 7 cr level in both fields. In this case the es- 
timate of the stellar mass depends only mildly on the assumed 
reddening, and therefore the correlation is likely to be intrin- 
sic, with more massive galaxies being also more absorbed. 

Given the previous two correlation, not surprisingly we 
also find a correlation between SFR and stellar mass (Fig- 
ure andll2F). The Spearman rank correlation coefficient 
is = 0.30 for Daddi-F, and r^ ~ 0.45 for Deep3a-F, indi- 
cating that the correlation between SFR and stellar mass is 
significant at the > 4<j level in both fields. The origin of the 
sharp edge in Fig. 112b ) and Fig. I12F ) is caused by the color 
limit BzK > —0.2 for selecting sBzKs. To show this clearly, 
we plot galaxies with —0.2 < BzK < 0.0 as open squares in 
Fig.ll2F). However, no or very few z ~ 2 galaxies exist below 
the BzK > —0.2 line, and therefore the upper edge shown in 
the figure appears to be intrinsic, showing a limit on the max- 
imum SFR that is likely to be present in a galaxy of a given 
mass. 

At z = the vast majority of massive galaxies 
(M* ^ IO^'Mq) are passively evolving, "red" galaxies (e.g., 
Baldry et al. 2004), while instead at z ~ 2, actively star- 
forming (sBzKs) and passive (pBzKs) galaxies exist in similar 
numbers, and Fig.^jshows that the most massive sBzKs tend 
also to be the most actively star forming. This can be seen 
as yet another manifestation of the downsizing effect (e.g., 
Cowie et al. 1996; Kodama et al. 2004; Treu et al. 2005), with 
massive galaxies completing their star formation at an earlier 
epoch compared to less massive galaxies, which instead have 
more prolonged star formation histories. 

Because of the correlations discussed above, UV-selected 
samples of z 2 galaxies (Steidel et al. 2004) will tend to 
preferentially miss the most star-forming and most massive 
galaxies. Still, because of the large scattering in the correla- 
tions, some of the latter galaxies will also be selected in the 
UV, as can be seen in Fig. 1121 and as emphasized in Shapley 
et al. (2004; 2005). 

6.4. Mass and SFR densities 

In this subsection we derive the contribution of the sBzKs 
to the integrated star formation rate density (SFRD) at z ^ 2 
and of sBzKs and pBzKs to the stellar mass density at z ^ 2. 
Some fraction of the sBzKs galaxies are known to be AGN- 
dominated galaxies (Daddi et al. 2004b; 2005b; Reddy et al. 
2005). To estimate the AGN contamination, we have used the 
80 ks XMM-Newton X-ray data that are available for Daddi- 
F (Brusa et al. 2005). A circular region of 11' radius from 
the point of maximum exposure time include 70 sBzKs 18 of 
which are identified with X-ray sources using a 5" radius error 
circle (see Fig. I12> . This fraction is comparable to the one 
estimated in the CDFS field (Daddi et al. 2004b) and in the 
GOODS-N field (Daddi et al. 2005b; Reddy et al. 2005), for 
which > 1 Ms Chandra data are available. Based also on the 
latter result, we assume the AGN contamination is about 25%, 
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Fig. 1 1 . — Reddening, star formation rate, and stellar mass histogram of sBzKs in Daddi-F and Deep3a-F. (a), (c), (e) (left panels): Plots for Daddi-F; (b), (d), 
and (f) (right): Plots for Deep3a-F. The filled area is the histogram for sBzKs in Daddi-F, which are associated with X-ray sources (about 25%). The dashed lines 
in (e) and (f) are the stellar mass histograms of pBzKs. [See the electronic edition of the Journal for the color version of this figure.] 



and we adopt this fraction to statistically correct properties 
derived from our sBzKs samples. 

The left panel of Fig. shows the differential contribu- 
tion to the SFR density at z 2± 2 from sBzKs as a function 
of their A'-band magnitude. Using the volume in the redshift 
range 1.4 ;S z ^ 2.5 (Daddi et al. 2004a; see also Reddy et 
al. 2005), an SFRD of 0.08 yr"' Mpc^^ is derived from 
the sBzKs (/Tvega < 20) in Deep3a-F, and an SFRD of 0.024 
Mq yr^' Mpc^^ is derived from the sBzKs (AVega < 19.2) in 
Daddi-F. These estimates are reduced, respectively, to 0.06 
Mq yr-1 Mpc-3 (/Tvega < 20) and 0.018 Mq yr"' Mpc^^ 
(^Tvega < 19.2) when subtracting an estimated 25% AGN con- 
tamination. Using the same method and the 24 sBzKs in 
the K20/GOODS-S sample, Daddi et al. (2004a) derived an 
SFRD 0.044 ± 0.008 Mq yr"' Mpc"^^ for the volume in the 
redshift range 1.4 <, z 2.5; i.e., ^ 25% lower than that de- 
rived here, possibly due to cosmic variance. However, note 
that there appears to be just a hint for the increasing trend in 
SFRD with K magnitude to flatten out at /Tvega ^ 20, indicat- 
ing that a substantial contribution to the total SFRD is likely 
to come from /Tvega > 20 sBzKs, and therefore the values de- 
rived here should be regarded as lower limits. 

Recently, Reddy et al. (2005) provided an estimate for the 
SFRD of TiTvega < 20 sBzKs of < 0.02 M© yr"' Mpc^^ about 



one-third of the estimate that we have obtained here. Part of 
the difference is due to the absence of the reddest sBzKs from 
the Reddy et al. (2005) sample, as already noticed. However, 
most of the difference is likely due to the fact that the Reddy 
et al. (2005) SFR estimate is based primarily on the X-ray 
emission interpreted with the Ranalli et al. (2003) relation. 
As shown by Daddi et al. (2005b), the X-ray emission inter- 
preted in this way typically underestimates the SFR of sBzKs 
by factors of 2-3, with respect to the radio-, mid-lR- and far- 
IR-based SFR estimates, all of which are also in reasonable 
agreement with the UV-corrected SFR estimate. 

The right panel of Fig. shows the differential contribu- 
tion to the stellar mass density at z ~ 2 from sBzKs and 
pBzKs as a function of their /T-band magnitude. The open 
squares and triangles represent values that were calculated 
from all sBzKs, the solid symbols represent values corrected 
for the AGN contamination. The stars and crosses represent 
the mass density contributed by pBzKs. 

The stellar mass density in Deep3a-F, integrated to our 
TiTvega < 20 catalog limit, is logp* =7.7 Mq, Mpc^^, in excel- 
lent agreement with the value reported in Table 4 of Fontana 
et al. (2004), i.e., logp* = 7.86 Mq Mpc"^ for 1.5 < z < 2.0 
galaxies, and logp* = 7.65 Mq Mpc^^ for 2.0 <z< 2.5 galax- 
ies, but now from a much bigger sample. These estimates 
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also agree with the logp* ^ 7.5 Mp, Mpc^^ estimate at z 2 
by Dickinson et al. (2003), although their selection is much 
deeper (Hab < 26.5), although it extends over a much smaller 
field (HDF). So, while our sample is likely to miss the con- 
tribution of low-mass galaxies, the Dickinson et al. sample is 
likely to underestimate the contribution of high-mass galaxies 
due to the small field and cosmic variance. There is little evi- 
dence for flattening of logp* by ^Tvega = 20. As already noted, 
the total stellar mass density at z ^ 2 has to be significantly 
larger than that estimated here, i.e., only from the contribu- 
tion of /Tvcga < 20 Bz/T-selected galaxies. 

There are 121 pBzKs in Deep3a-F, and for ^ 100 of 
them we derive > lO'^ Mg,. Correspondingly, the num- 
ber density of pBzKs with > 10^' Mq over the range 
1.4<z<2.0 is (1.8 ±0.2) x lO"'* Mpc"'^ (Poisson error 
only). This compares to 3.4 x 10^^ Mpc^-' over the same 
redshift range as estimated by Daddi et al. (2005a) using six 
objects in the Hubble Ultra Deep Field (HUDF) with spectro- 
scopic redshift. While the Daddi et al. (2005a) HUDF sample 
is important to establish that most pBzKs are indeed passively 
evolving galaxies at 1 .4 < z < 2.5, their density measurements 
is fairly uncertain due to cosmic variance. Being derived from 
an area which is ~ 30 times larger than HUDF, the results 
presented here for the number density of massive, passively 
evolving galaxies in Deep3a-F in the quoted redshift range 
should be much less prone to cosmic variance. Hence, we 
estimate that, compared to the local value at z (9 x lO^'* 
Mpc"^, Baldry etal. 2004), at 1.4 < z < 2 there appears to 
be about 20%±7% of massive (> lO'^M©), passively evolv- 
ing galaxies, with the error above accounting also for cosmic 
variance. 

7. SUMMARY AND CONCLUSIONS 

This paper presents the results of a survey based on BRIzJK 
photometry obtained by combining Subaru optical and ESO 
near-IR data over two separate fields (Deep3a-F and Daddi- 
F). Complete K-selected samples of galaxies were selected 
to Kycga < 20 in the Deep3a-F over 320 arcmin-^, and to 
^Vcea ~ 19 in the Daddi-F over a field roughly twice the area. 
Deep multicolor photometry in the BRIz bands were obtained 
for the objects in both fields. Object catalogs constructed 
from these deep data contain more than lO'* objects in the NIR 
bandpasses. Galaxy A'-band number counts were derived and 
found to be in excellent agreement with previous survey re- 
sults. 

We have used color criteria to select candidate mas- 
sive galaxies at high redshift, such as Bz^T-selected star- 
forming (sBzKs) and passively evolving (pBzKs) galaxies at 
1.4 ^ z 2.5, and EROs, and derived their number counts. 
The main results can be summarized as follows. 

1. Down to the /T-band limit of the survey the log of the 
number counts of sBzKs increases linearly with the K magni- 
tude, while that of pBzKs flattens out by Ky^ga ~ 19. Over 
the Deep3a-F we select 387 sBzKs and 121 pBzKs down 
to ^^vega = 20, roughly a factor of 10 more than over the 
52 arcmin^ fields of the K20 survey. This corresponds to 
a ^ 30% higher surface density, quite possibly the result of 
cosmic variance. Over Daddi-F we select 108 sBzKs and 48 
pBzKs down to /Tvega = 19.2. 

2. The clustering properties (angular two-point correlation 
function) of EROs and Bz^T- selected galaxies (both sBzKs and 
pBzKs) are very similar, and their clustering amplitudes are 
about a factor of 10 higher than those of generic galaxies in 
the same magnitude range. The most strongly clustered pop- 



ulations at each redshift are likely to be connected to each 
other in evolutionary terms, and therefore the strong cluster- 
ing of EROs and BzKs makes quite plausible an evolutionary 
link between BzKs at z ~ 2 and EROs at z ^ 1, with star for- 
mation in sBzKs subsiding by z ~ 1 thus producing passively 
evolving EROs. While some pBzKs may well experience sec- 
ondary, stochastic starbursts at lower redshift, the global evo- 
lutionary trend of the galaxy population is dominated by star 
formation being progressively quenched in massive galaxies, 
with the quenching epoch of galaxies depending on environ- 
mental density, being earlier in high-density regions. 

3. Using approximate relations from Daddi et al. (2004a) 
and multicolor photometry, we estimated the color excess, 
SFR and stellar mass of sBzKs. These /Tvega < 20 giiliix- 
ies have median reddening E{B — V) ^ 0.44, average SFR^ 
190 M0>T^', and typical stellar masses ^ lO^'Mr:.. Correla- 
tions between physical quantities are detected: the most mas- 
sive galaxies are those with the largest SFRs and optical red- 
dening E{B — V). The high SFRs and masses of these galaxies 
add further support to the notion that these z — 2 star-forming 
galaxies are among the precursors of z — 1 passive EROs and 
z ~ early-type galaxies. 

4. The contribution to the total star formation rate density 
at z ~ 2 was estimated for the AVega < 20 sBzKs in our fields. 
These vigorous starbursts produce an SFRD ^ 0.06 Mq yr^' 
Mpc^^, which is already comparable to the global SFRD at 
z ~ 2 as estimated from other surveys and simulations (e.g. 
Springel & Hernqwist 2003; Heavens et al. 2004). However, 
a sizable additional contribution is expected from AVega > 20 
sBzKs. 

5. In a similar fashion, the stellar mass of pBzKs was ob- 
tained, with the result that the number density of /Cvega < 20 
pBzKs more massive flian IO^^Mq is about 20%±7% of that 
of similarly massive, early-type galaxies at z = 0, indicat- 
ing that additional activity and subsequent quenching of star- 
formation in >, 10^ 'M© star-forming galaxies must account 
for increasing the number of massive passive galaxies by 
a factor of about 5 from z = 1.7. The number density of 
^ 10^ sBzKs is similar to that of pBzKs. Given flieir 
strong star-formation activity, it seems that by z ^ 1-1.4 the 
full population of local ^ 10^ 'M© passive galaxies could be 
eventually assembled as a result. 

This result, advocated also in Daddi et al. (2005b), may 
appear in contradiction with the recent finding by Bell et al. 
(2004) of a factor of 2 decrease in the number density of early- 
type galaxies at z ~ 1, with respect to the local value (see 
also Faber et al. 2005). However, our analysis of the BeU 
et al (2004) results shows that most of this evolution is to be 
ascribed to the progressive disappearence with increasing red- 
shift of the fainter galaxies, while the population of the bright- 
est, most massive galaxies being substantially stable. This 
would be, in fact, another manifestation of the downsizing ef- 
fect. A future publication will address this point in its full 
details and implications. 

Mapping the metamorphosis of active star-forming galax- 
ies into passively evolving, early-type galaxies from high to 
low redshifts, and as a function of galaxy mass and environ- 
ment is one of the primary goals of the main ongoing galaxy 
surveys. Using Subaru and VLT telescopes, optical and near- 
infrared spectra are being obtained, with targets from the 
present database having been selected according to the same 
criteria adopted in this paper. Future papers in this series will 
present further scientific results from this pilot survey, along 
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with a variety of data products. 

We thank the anonymous referee for useful and construc- 
tive comments that resulted in a significant improvement of 
this paper The work is partly supported by a Grant-in- Aid 
for Scientific Research (16540223) by the Japanese Ministry 
of Education, Culture, Sports, Science and Technology and 



the Chinese National Science Foundation (10573014). X.K. 
gratefully acknowledges financial support from the JSPS. 
E.D. acknowledge support from NASA through the Spitzer 
Fellowship Program, under award 1268429. L.F.O. thanks 
the Poincare fellowship program at Observatoire de la Cote 
d' Azur and the Danish Natural Science Research Council for 
financial support. 



REFERENCES 



Abraham, R. G., et al. 2004, AJ, 127, 2455 

Adelberger, K. L., & Steidel, C. C. 2000, ApJ, 544, 218 

Adelberger, K. L., Erb, D. K., Steidel, C. C, et al. 2005, ApJ, 620, L75 

Baldry, I. K., Glazebrook, K., Brinkmann, J., et al. 2004, ApJ, 600, 681 

Bell, E. E, et al. 2004, ApJ, 608, 752 

Bertin, E. & Amouts, S. 1996, A&AS, 117, 393 

Broadhurst, T. J., Ellis, R. S., & Glazebrook, K. 1992, Nature, 335, 55 

Brown, M. J. I., Jannuzi, B. T., Dey, A., Tiede, G. R 2005, ApJ, 621, 41 

Brusa, M., et al. 2005, A&A, 432, 69 

Bmzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Calzetti, D., 2004, proceedings of the conference 'Starbursts - from 30 
Doradus to Lyman Break galaxies', Cambridge, Sept 2004, eds. de Grijs 
and Gonzalez-Delgado i astro-ph/0501027 i 

Calzetti, D., Armus, L., Bonun, K. i^., et al. 2000, ApJ, 533, 682 

Cardelli, J. A., Clayton, G. C, Mathis, J. S. 1989, ApJ, 345, 245 

Carlberg, R. G., Cowie, L. L., Songaila, A., Hu, E. M. 1997, ApJ, 484, 538 

Cimatti, A., et al. 2002a, A&A, 392, 395 

Cimatti A., et al. 2002,A&A, 381, L68 

Cimatti, A., et al. 2004, Nature, 430, 184 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, MNRAS, 319, 
168 

Cowie, L.L., Songaila, A., Hu. E.M., & Cohen, J.G. 1996, AJ, 112, 836 
Cutri, R. M., et al. 2003, VizieR Online Data Catalog, 2246, 
Daddi, E., Cimatd, A., Pozzetti, L., et al. 2000, A&A, 361, 535 
Daddi, E., Broadhurst, T., Zamorani, G., et al. 2001, A&A, 376, 825 
Daddi, E., et al. 2003, ApJ, 588, 50 

Daddi, E., Cimatd, A., Renzini, A., et al. 2004a, ApJ, 617, 746 
Daddi, E., et al. 2004b, ApJ, 600, L127 
Daddi, E., Renzini A., Pirzkal N., et al. 2005a, ApJ, 626, 680 
Daddi, E., Dickinson, M., Chary, R., et al. 2005b, ApJ, 631, L13 
Dickinson, M., Papovich, C, Ferguson, H. C, & Budavari, T. 2003, ApJ, 
587, 25 

Doherty, M., Bunker, A. J., Ellis, R. S., McCarthy, P J. 2005, MNRAS, 361, 
525 

Elston, R., Rieke, G. H., & Rieke, M. J. 1988, ApJ, 331, L77 
Faber, S. M., et al., 2005, submitted to ApJ i astro-ph/0506044l 
Fang, E, et al. 2004, ApJS, 154, 35 

Ferguson, H. C, Dickinson, M., & WilHams, R. 2000, ARA&A, 38, 667 
IFirth, A. E., et al. 2002, MNRAS, 332, 617 
Fontana, A., et al. 2004, A&A, 424, 23 
Franx, M., et al. 2003, ApJ, 587, L79 

Gardner, J. R, Cowie, L. L., & Wainscoat, R. J. 1993, ApJ, 415, L9 
Georgakakis, A., Afonso, J., Hopkins, A. M., et al. 2005, ApJ, 620, 584 
GiavaHsco, M., et al. 2004, ApJ, 600, L93 

Glazebrook, K., Ellis, R., Santiago, B., & Griffiths, R. 1995, MNRAS, 275, 
L19 

Glazebrook, K., et al. 2004, Nature, 430, 181 

Granato, G. L., De Zotti, G., Silva, L., et al. 2004, ApJ, 600, 580 

Heavens, A., Panter, B., Jimenez, R., Dunlop J. 2004, Nature, 428, 625 

Huang, J.-S., et al. 2001, A&A, 368, 787 

Kauffmann, G. & Chariot, S. 1998, MNRAS, 297, L23 



Kerscher, M., Szapudi, I., & Szalay, A. S. 2000, ApJ, 535, L13 
Kodama, T., Yamada, T., Masayuki, A., et al. 2004, MNRAS, 350, 1005 
Kong, X., Chariot, S., Brinchmann, J., & Fall, S. M. 2004, MNRAS, 349, 
769 

Kummel, M. W. & Wagner, S. J. 2001, A&A, 370, 384 

Landolt, A. U. 1992, AJ, 104, 372 

Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 

Lejeune, T., Cuisinier, E, & Buser, R. 1997, A&AS, 125, 229 

McCarthy, R J. et al. 2001, ApJ, 560, LI 31 

McCarthy, R J. 2004, ARA&A, 42, 477 

McCarthy, R J., et al. 2004, ApJ, 614, L9 

McCracken, H. J., Metcalfe, N., Shanks, T., et al. 2000, MNRAS, 311, 707 
McLeod, B. A., Bernstein, G. M., Rieke, M. J., et al. 1995, ApJS, 96, 1 17 
Menci, N., Cavaliere, A., Fontana, A., et al. 2002, ApJ, 575, 18 
Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64 
Miyazaki, S. et al. 2002, PASJ, 54, 833 

Moorwood, A., Cuby, J.G., Lidman, C. 1998, The ESQ Messenger 91,9 
Moustakas, L. A., Davis, M., Graham, J. R., et al. 1997, ApJ, 475, 445 
Moustakas, L. A., et al. 2004, ApJ, 600, L131 

Nagamine, K., Fukugita, M., Cen, R., et al. 2001, MNRAS, 327, LIO 
Nagamine, K., Cen, R., Hernquist, L., et al. 2005, ApJ, 627, 608 
Oliver, S., et al. 2004, ApJS, 154, 30 
Ouchi, M., et al. 2004, ApJ, 61 1, 660 

Persson, S. E., Murphy, D. C, Krzeminski, W., et al. 1998, AJ, 116, 2475 
Pickles, A. J. 1998, PASP, 110, 863 

Ranalli, P, Comastri, A., Setti, G. 2003, A&A, 399, 39 
Reddy, N. A., Erb, D. K., Steidel, C. C, et al. 2005, ApJ, 633, 748 
2Renzini, A., & da Costa, L. N. 1997, The Messenger, 87, 23 
Roche, N., Eales, S. A., Hippelein, H., et al. 1999, MNRAS, 306, 538 
Roche, N. D., Dunlop, J., & Almaini, O. 2003, MNRAS, 346, 803 
Saracco, P, lovino. A., Garilh, B., et al. 1997, AJ, 1 14, 887 
Saracco, P, et al. 2005, MNRAS, 357, L40 

Schlegel, D. J., Finkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 
Shapley, A. E., Erb, D. K., Pettini, M., et al. 2004, ApJ, 612, 108 
Shapley, A. E., Steidel, C. C, Erb, D. K., et al. 2005, ApJ, 626, 698 
Smith, J. A., et al. 2002, AJ, 123, 2121 

Somerville, R. S. 2004a, in Multiwavelength Mapping of Galaxy Formation 

and Evolution, eds. R. Bender and A. Renzini, astro-ph/0401570 
Somerville, R. S., Lee, K.. Ferguson, H. C, et al. ZUU4B, Apj, ouu, L171 
Somerville, R. S., et al. 2004c, ApJ, 600, L135 
Springel, V., Hernquist, L. 2003, MNRAS, 339, 312 
Steidel, C. C, Giavalisco, M., Pettini, et al. 1996, ApJ, 462, L17 
Steidel, C. C, Shapley, A. E., Pettini, M., et al. 2004. ApJ, 604, 534 
Thompson, D., et al. 1999, ApJ, 523, 100 

Treu, T, Ellis, R.S., Liao, T.X., & Cohen, J.G. 2005, ApJ, 622, L5 
Vandame, B. 2002, Proc. SPIE, 4847, 123 

Yagi, M., Kashikawa, N., Sekiguchi, M., et al. 2002, AJ, 123, 66 
Yan, L., Thompson, D., Soifer, B. T. 2004, AJ, 127, 1274 



